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Abstract:  Models  for  non-linear  dielectrics  and  magnetic  ferrites  are  developed  and  coded  into 
the  particle-in-cell  code  ICEPIC.  The  non-linear  dielectric  model  includes  a  relaxation  term  to 
account  for  the  physical  response  time  of  the  material,  which  may  be  longer  than  an  FDTD  time 
step.  The  ferrite  model  accounts  for  the  non-linearity  of  the  Landau-Lifshitz-Gibert  equation, 
and  the  magnetization  update  is  derived  from  the  solution  of  the  equation  on  a  time  step.  Initial 
testing  of  the  models  is  encouraging,  and  work  is  underway  on  further  improvement  and  refinement. 

Keywords:  FDTD,  non-linear 


1.  Introduction 

The  improved  concurrent  electromagnetic  particle-in-cell  (ICEPIC)  code  implements  finite-difference 
time-domain  (FDTD)  electromagnetics  and  the  relativistic  Lorentz  force  law  for  tracking  charged 
particles.  The  code  is  massively  parallel  and  applicable  to  modeling  a  variety  of  high  power  mi¬ 
crowave  (HPM)  devices  [1],  [2],  More  recently,  there  is  interest  in  the  use  of  non-linear  materials  in 
the  generation  of  HPM.  Two  types  of  materials  are  of  interest,  non-linear  dielectric  materials  and 
ferrite  materials,  which  are  governed  by  the  Landau-Lifshitz-Gilbert  equation.  Although  previous 
authors  consider  similar  materials  [3]- [6],  to  our  knowledge,  aspects  of  our  formulation  and  concerns 
related  to  the  generation  of  high  power  microwaves  are  unique. 


2.  Formulation 


A.  Non-linear  Dielectrics 

Non-linear  dielectrics  are  materials  for  which  the  permittivity  depends  on  the  electric  field  strength. 
Of  particular  interest  are  materials  for  which  the  permittivity  is  a  decreasing  function  of  the  elec¬ 
tric  field.  In  such  materials,  higher  field  values  experience  lower  permittivity,  hence  faster  wave 
propagation  speed.  This  leads  to  a  steepening  of  a  wave  pulse  as  it  propagates,  or  in  other  words, 
the  formation  of  a  shock  wave. 

The  approach  to  modeling  non-linear  dielectric  material  is  to  curve  fit  an  analytical  function 
to  experimental  permittivity  data.  Two  fitting  functions  considered  are  the  Lorentzian  function, 
given  by 


A  o.5r 

^  ’  ~  7T  (E  -  E0)2  +  (0.5T)2 


(1) 
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(a)  Barium  Strontium  Titanate  (b)  Lead  Magnesium  Niobate  38  (c)  Lead  Magnesium  Niobate  85 


Figure  1:  Experimental  relative  permittivity  and  associated  curve  fits  for  materials  of  interest. 


where  A,  Eq,  and  T  are  parameters  of  the  fit.  The  primary  disadvantage  to  the  Lorentzian  function 
is  er(E  — >  oo)  — *  0;  permittivity  values  below  one  are  problematic.  Thus,  the  second  (sixth  order) 
fitting  function,  given  by 


tr(E)  —  £r,o o  T 


er  0  ^r 


1  +  a\E2  +  a2E4  +  CI3E6 


(2) 


where  eri0o,  er, o,  at,  02,  and  03  are  parameters  of  the  fit.  For  this  function,  er(0)  =  er,o  and 
er(E  —>  00)  er)00.  Plots  of  experimental  data  and  associated  curve  fits  for  three  materials  of 
interest  are  shown  in  Fig.  1.  Note  that  in  addition  to  being  non-linear,  the  permittivity  values  in 
these  materials  are  quite  large. 

The  incorporation  of  the  non-linear  permittivity  into  Maxwell’s  equations  involves  the  inter¬ 
polation  of  the  staggered  electric  field  values  from  the  previous  time  step  to  find  the  magnitude 
of  E  at  a  given  held  point,  then  the  computation  of  the  permittivity  based  on  E.  The  computed 
permittivity  is  then  used  to  update  first  D,  then  E  from  Ampere’s  law  and  from  D  =  eE.  For 
stability,  it  is  important  that  only  past  values  of  E  contribute  to  the  computed  permittivity  [7]. 

The  physical  response  time  of  a  material  is  often  much  slower  than  an  FDTD  time  step.  Thus, 
in  the  formulation,  the  response  of  the  material  is  incorporated  using 

a^  +  €r  =  f{E)  (3) 


where  f(E)  is  er(E)  from  Eq.  1  or  2.  Analytically  solving  this  equation  on  a  time  step  yields 


e?  =  e"-1e~aAt  + 


jjjn—1  _  j^jn—2 


E71-1 


[1  -  e~aAt] 


df 

+  Af  — 
*  8E 


j^n—1  En~ 2 


E71-1 


(4) 


B.  Magnetic  Ferrites 
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In  many  cases,  a  ferrite  can  be  modeled  by  assuming  a  large,  static  magnetic  field  and  forming 
a  linearized  gyrotropic  permeability  [8].  However,  the  applications  of  interest  do  not  admit  such  a 
model.  Thus,  the  magnetization  is  modeled  using  the  Landau-Lifshitz-Gilbert  equation, 


<9M 

~dt 


7o 

1  +  A2 


x  Heff  + 


A 

M 


M  x  M  x  H 


where  7o  =  2.212761045  x  105m/(A  s), 


(5) 


Hcff  =  H  +  H, 


2  kn 


static 


/io|M|: 


(M  -u)u  + 


2  A 


//q|M| 


rV2M 


(6) 


B  =  M  (H  +  M)  (7) 

and  |M|  =  Ms,  A  (sometimes  denoted  a),  ku,  u,  and  A  are  parameters  of  the  material  being 
modeled.  Note  that  Eq.  5  keeps  the  magnitude  of  M  constant  with  time  but  varies  its  direction. 

Ferrite  materials  are  modeled  by  placing  all  components  of  the  magnetization  vector  (M)  at 
the  center  of  a  FDTD  cell.  The  magnetic  field  values  are  interpolated  from  the  staggered  Yee  grid 
to  the  center  of  the  cell,  and  the  quantity  Heff  =  Heg-  +  (A/|M|)M  x  Heff  is  defined.  Eq.  5  then 
becomes 


DM 

~dt 


7o 

1  +  A2 


M  X  Heff 


(8) 


which  can  be  solved  on  t  =  [{n  —  0.5)Af,  (n  +  0.5)Aj)  to  give 
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(9) 


In  practice,  is  approximated  with  B  is  updated  using  Faraday’s  law,  then  the  M 

values  are  interpolated  to  the  H  locations,  and  H  is  updated  from  Eq.  7. 

3.  Results 

The  dielectric  model  is  tested  by  simulating  a  parallel  plate  capacitor  filled  with  non-linear 
dielectric.  A  known  current  is  driven  from  one  plate  to  the  other,  and  the  electric  held  between  the 
plates  is  observed.  From  Gauss’  law,  if  A  is  the  area  of  the  plates,  respectively,  the  current  from 
one  plate  to  the  other  is  given  by 

j  -  di = 4 |£(E)E1  (10) 
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(a)  Barium-Strontium-Titanate  (b)  Lead-Magnesium-Niobate  38  (c)  Lead-Magnesium-Niobate  85 


Figure  2:  Comparison  of  computed  current  to  drive  current  for  parallel  plate  capacitor  filled  with 
three  different  non-linear  dielectrics.  For  this  test,  the  plate  area  is  2.5  mm2,  the  plate  separation 
is  0.632mm,  and  a  =  lOps. 


(a)  Parallel  propagation,  Ihc-pol  (b)  Parallel  propagation,  rhc-pol 


Frequency  (GHz) 


Frequency  (GHz) 


(c)  Perpendicular  propagation,  (d)  Perpendicular  propagation, 
parallel  magnetic-pol  perpendicular  magnetic-pol 


Figure  3:  Propagation  characteristics  for  a  er  =  10,  a  =  0,  Ms  =  159kA/m,  A  =  0.05,  ku  =  0, 
A  =  0  ferrite  in  a  71.0kA/m  static  magnetic  field. 


The  current  computed  from  Eq.  10  is  compared  to  the  drive  current.  Results  are  shown  in  Fig.  2, 
and  the  agreement  is  good. 

The  ferrite  model  is  tested  by  evaluating  wave  propagation  in  a  ferrite  in  a  large  static  magnetic 
held.  In  this  regime,  the  ferrite  behaves  as  a  linearized  gyrotropic  material,  and  the  wave  number 
can  be  analytically  computed  [8].  Analytic  and  simulated  wave  number  versus  frequency  plots  are 
shown  in  Fig.  3  for  a  ferrite  with  er  =  10,  a  =  0,  Ms  =  159kA/m,  A  =  0.05,  ku  =  0,  A  =  0  in  a 
71.0kA/m  static  magnetic  held.  Results  are  shown  for  left-hand  and  right-hand  circular  polarization 
when  the  propagation  is  parallel  to  the  magnetic  held,  and  for  the  magnetic  held  polarized  parallel 
and  perpendicular  to  the  static  magnetic  held  when  the  propagation  is  perpendicular  to  the  static 
held.  Except  when  the  attenuation  of  the  ferrite  is  large  (2-8  GHz  for  the  right-hand  circular 
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polarized  case  and  4-8  GHz  for  the  perpendicular  polarized  case),  the  agreement  is  excellent. 

In  both  models,  a  slow  growth  mode  is  sometimes  observed.  Investigation  shows  that  the 
growth  is  related  to  aliasing  when  the  wave  number  becomes  larger  than  the  Nyquist  cut-off  of  the 
numerical  grid.  Work  is  underway  to  improve  the  utility  of  the  models  by  suppressing  this  mode 
[9].  Results  are  to  be  reported  in  the  conference  presentation. 


4.  Conclusion 

Models  for  non-linear  dielectrics  and  magnetic  ferrites  are  developed  and  coded  into  the  particle- 
in-cell  code  ICEPIC.  The  non-linear  dielectric  model  includes  a  relaxation  term  to  account  for  the 
physical  response  time  of  the  material,  which  may  be  longer  than  an  FDTD  time  step.  The  ferrite 
model  accounts  for  the  non-linearity  of  the  Landau-Lifshitz-Gibert  equation,  and  the  magnetization 
updated  is  derived  from  the  solution  of  the  equation  on  a  time  step.  Initial  testing  of  the  models 
is  encouraging,  and  work  is  underway  on  further  improvement  and  refinement. 
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